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I. INTRODUCTION 


During the past few decades, weather forecasting accuracy has increased largely due 
to the improvements in numerical modeling of the atmosphere. Forecasts of 1 - 3 days can 
be made with a good deal of certainty, allowing for necessary preparations to be made for 
oceanic or continental winter storms. However, rapidly changing events such as explosive 
cyclogenesis are still often missed by the numerical models and the human forecaster. 
This remains a potentially hazardous situation since people are often caught unprepared 
with little or no warning. 

While the fundamentals of cyclogenesis are generally understood and documented, the 
factors which combine to create a rapidly deepening cyclone are still unclear. Opertitional 
forecast models used for 1 to 3 d:iy forecasts have consi.stently tended to under-forecast the 
deepening rate of these storms (Sanders and Gyakum 1980), although more recent studies 
have shown an improvement in this area (Sanders 1987). A number of reasons have been 
suggested to account for this failure including spatial resolution, inititd conditions, and 
subgrid-scale parameterizations. 

.A large number of studies using numerical model simulations of explosive cyclogene¬ 
sis have been conducted. In general, these studies fall into two categories. The first group 
includes near perfect model simuUitions allow'ing for diagnostic studies. These studies 
allow researchers to get a more complete and dynamically consistent picture ol the de¬ 
veloping cyclone, even though actual data measurements are sparse temporally and spatial¬ 
ly (e.g.. Whitaker et al. 1988; Uccellini et al. 1987). The second kind ol study involves 
investigations of model failures to match the ob.served deepening rate of an extratropical 
cyclone (e.g.. Anthes et al. 1983). As numerical models usually under-predict the ob¬ 
.served deepening rate, various parameterizations and resolutions are varied in an attemjii 
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to improve the model forecast, witii varying degree.-, of success. This thesis, however, will 
investigate a failure of a numerical model w'hich forecast overdevelopment of an explo¬ 
sively deepening extratropical cyclone. 

Recently, large scale experiments have been conducted in order to study the explosive 
deepening of extratropical cyclones. The most recent of the.se was the Experiment on 
Rapidly Intensifying Cyclones over the Atlantic (ERICA), conducted from 1 December 
1988 to 1 March 1989 (Hadlock 1988; Hartnett et al. 1989). Its purpose was to provide a 
large data sampling of .such cyclones over the open waters of the western North Atlantic 
Ocean. Aircraft, buoys, increased rawinsondes, and radar data supplemented the usual 
satellite imagery during Intense Observation Periods (lOP's). 

The second storm of the experiment, IOP'2, devek^ped on l.>-14 December 1988. The 
Naval Operational Regional Atmospheric Prediction System (NORAPS) forecast initial¬ 
ized ill 12Z on 12 Dec 88 missed the time and rate (d' rapid deepening of the IOP-2 cy¬ 
clone, although its 48 h forecast of centrtil pressure was only 6 mb off. Using additional 
buoy and drop.scnde data obtained during ERICA, this study will investigate the successful 
and unsuccessful as[)ects of the model's forecasts. 

A literature review' is presented in section 2 focusing on the hitent heating aspect of 
rtipid cyclogenesis. Section 3 summarizes the synoptic events of lOP-2 and briefly de¬ 
scribes the operational numerical forectists (d this cyclone. A brief description of the 
NORAPS model is given in section 4. The numcrictil simubitions mtide using the 
NORAPS model are documented in section followed by a diagnostic study of these 
simulations in section 6. 
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11. BACKGROUND 


Numerous studies have been conducted to investigate the causes of explosive cyclo¬ 
genesis and the reasons for poor forecasts of these events by numerical models. Anthes et 
al. (1983) described model simulations of the QE II storm. Their study concluded that 
lower tropospheric baroclinic instability was the major cause of growth for that cyclone. 
This result was confirmed by composite rawinsonde data studies from weather ships 
(Rogers and Bosart 1986). However, Uccellini (1986) has show n that upper-tropospheric 
conditions w'ere also important in the OE II storm. 

Anthes et al. (1983) state that latent heat release appeared to be mi^re important in the 
latter stages of development. Numerical simulations with and without latent heating 
show'ed little difference in the 12 h forecast, indicating the lack of importance in latent 
heating for the early period of development of that case. While the incretised latent betning 
improved the intensity forectist. the predicted motion of the sti'rm w tis slow ed due to this 
heating, resuning lu iaiger posiiniii eiioi.>. Recem .siudic:^ al.so li.ivc .>hov. n that the verlictil 
distribution of latent heating has a significtmt effect on the deepening the storm, w ith ;i 
low er height of mti.xim.um Iteating producing a greater intensification. Danard (16S5) 
postulates that this is due to a lower level of non-uiveigeri^c aiid ii.C!et...od ;..v. le'.e! 
convergence and vorticity production. The heating of the lower troposphere also causes 
destabilization, resulting in ;i more rtipid intensification. 

Danard (1985) studied two explosive deepening events in the eastern North Pacilic 
using an 8-level large scale (190 km grid spacing) primitive equation model. He too 
crmcluded that both latent and sensible heating play a major role in cyclone intensification. 
Three variations of the cumulus parameterization scheme developed by Kuo (1 074) 





indicated that the maximum amount of convective precipitation was directly propt'rtional 
to the amount of deepening, although the differences w'ere not large (5 mh at most). 

Leslie et al. (19iS7) conducted model simulations of an explosively deepening cyclone 
off the east coast of Australia. Their findings on the role of precipitation agreed with 
Anthes et al. (198.^). For their ca.se study it appears that cumulus convection is not impi'r- 
tant during the early sUiges of the cyclone formtition, but plays a decisive role in the rapid 
intensification stage. They also tested the importance (tf other small scale processes such 
as effects of topographv, surface fluxes, and sea-surface temperature analyses. The\ 
concluded that all effects \'were necessary for evclone development anti that none could be 
timitted without degrading the lorecast sigmfictmtly. Cven with 2.^ km hori/iuital lesolu- 
tion, they ftiiled t(' achie\e the t'bserved deepening rate. 

Mullen and Baumhefner (1988) pointed out th;it case studies <9 intii\ idual storms can 
be misleading ;md thtit the results ctinnot always be genertilized. 'I'heir studs used an 1 1 
western Ntartn Pacitic Ocetin storm ensemble frtmi the NCAR Community Climate .Motiel. 
From ;i L^d-day simuhition of perpetual Janutiry conditions, their results intlicate that total 
diabatic heating and baroclinic dynttmics contribute equtilly to storm deselopment. .Addi- 
tiontilly. sensible and Latent heating are etich responsible for approximtitels one-halt I't the 
deepening due to diabatic processes. They also found thtit the deepening rate v.is insensi¬ 
tive to the exact specdictition of precipiititio.i pt.rtin.eterization although the_\ camioned 
nt.. to rule out other sclicme.s. 

Robertson and Smith (198,’) conducted a diagnostic energetics tintilysis of extraiiopi- 
cal cyclone development using the Limited Area .Mesosctile Prediction Ssstem (L.WIPS) 
from Drexel Ums ersity. Two versions of the model, mc'ist and dry. uere run .m i\s o case 
studies. Their results also show that the baroclinic aspects of development are greatly 
enhanced bv both convective :md stable latent heat reletise in the model. In the moist 
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l\)recr;sts. a substantial amount ot' eddy kinetie energy is generated due tct ageostrophic 
fltnv induced by diabatic heating. 

A comprehensive study on numerical prediction ot explosiv e cyclogenesis u as recent¬ 
ly conducted by Kuo and Ltw-Nam (IWO) 'I hey ran 24 h simulations ot 9 nine explosive 
cyclones over the western Atlantic ocean, varying the grid resolution, cumulus parameteri¬ 
zation. initial and boundary conditions, as well as surtace fluxes. Their results indicate 
that the most crucial model component for short range (0-24 h) prediction of rapid cyclo¬ 
genesis is the initial conditions, followed by the iKui/iuttal grid resolution, precipil.ition 
parameterizatitm. and lateral birundary coiuiiiions. w iili surtace fluxes plav ing a minor 
role. More interesting is their finding thal upnclu convection actually degraded the fore¬ 
cast due tit an incorrect location of heating, (irid scale precipitation (i.e.. non-convective 
or stable) associtited with the slantwise ascent in the warm front was crucitil for rapid 
development. Simulations using the Arak.iwa'Scluibert (10~4) scheme produced superior 
forecasts to tla.'se using Kuo (19~4). Hut this wa'' thought to be ctiused t\v the fact that tile 
Arakawa-Schubert scheme resulted in much lower convective precifutation amounts, thus 
allowing mi'ie ot tlie tivailable moiMure to !'e usevi hv tlie stable |irecipitation scheme. 
Thev also concluded that the Kuo scheme w,ix too ea-'V to aetivate. 

Most recentlv. twn studies w ere coiulucted I'H cvcioiu"' from the Clenesis ot .Aihintic 
Lows [-.xperiment ((l.Al.L) using the NOK AI’S model, l.iou et al. ( Ltooi ran the mod.el on 
the Iniensiv e Ot’serv alion PeriiHl (K )P) 2 c\ elone w ilh and w ithoiil tlie extra data obtained 
as ;i result of experiment. Since the eveiogenesis and dev i-loj'ment were along the east 
coast of the I’.S.. the additional dat.i trom the exiiermient resulted in onlv a limited im- 
jiiovement on the NOl^APS analyses and lorecasts le.g.. 1 - mb in >,entral pressure). In 
general, NOR.APS predicted the eyelogenesis event lairlv well in the short term (24 - .>(i h) 
but missed the strength tinii timing of a cold surue resulting in poiir forecasts later in the 
period. Seeontlary eyelogenesis at .’oi li was eorreeilv lorecast bv the model. Wash el al. 





(IWO) studied the GALE IOP9 cyclone. Using the additional GALE observations, the 
NORAPS optimal interpolation analysis reveals critical subsyni^ptic features which were 
important in the development. However, NORAPS incorrectly deepened the wrong low 
due to a poorly forecast 250 mb jet streak. Central pressures were under-forecast by 10-15 
mb after 36 h. NORAPS forecasts using the additional GALE data were better on the 
storm track but still deficient on intensity. 

The current study also investigates the importance of latent heating in the develop¬ 
ment of an explosively deepening cyclone. However, the regional model used in this study 
operationally over-predicted the deepening of the ERIC.A lOP-2 cyclone. The importance 
of the model initial conditions also will be studied. ERICA data (aircraft dropsondes. 
rawinsondes. and buoys) will be used to verify the model performance. The results should 
yield a better understtmding of numerical modeling processes in forecasts of explosive 
cyclone de\’elopment. 
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III. ERICA IOP-2 SYNOPTIC SUMMARY 


During the development of the ERICA lOP-2 cyclone, two disturbances of different 
origins interacted over the western North Atlantic Ocean to produce a rapidly deepening 
surface low (Chalfant 1989). Figure 1 presents the Fleet Numerical Oceanography Cen¬ 
ter’s (FNOC) 500 mb analysis for 1200 UTC 12 December 1988^. A large trough lies 
over eastern Canada and extends along the east coast (d' the United States. East of the 
Virginia coastline, two air streams merge, one origin:), ng from Northwest Canada, the 
other from the southeastern U.S. 

By 1.5/0000, short wave troughs are apparent in each td’ the air streams (Fig. 2 ). The 
northerly wave is located over Michig;in and Wisconsin with the southerly disturbance 
over the southetist U.S. Note the phase l:ig (d' the -.'^0"C isotherm compared to the height 
and wind fiel.i in the northern trough (o\er Minnes))t;i). indic;iting the likelihood of addi- 
tiontil barcK’linic development. The FNOC' surface an;i!\sis for this time (Fig. 2 ) reveals a 
cold front extending southwestward from a h'w near 40°N 50'W to the southern tip of 
Florida. Behind it. a high pressure system is nun ing cdi the New England coastline. The 
V. inds associated with the high are advecting cold air over the w arm Gulf Stream waters. 
Note that buoy 41001 ;it .55°N is reporting an air temjierature (d lO'C and a sea- 

surface (d' 20.5"C. 

A surface low begins to de\ elop ahnig the tail emi td the cold front off the Florida 
coast after l.i/OOOO in response to the southerly short wa\e. The low slowly deepens and 
moves northeastw ard over the next 12 h as seen in Fig. 4. Meanw hile at 500 mb. the 
northerly trough has continued to develop and move to the Ohio Valley (Fig. 5). Note. 

.. Ilereaf'ler. dale and lime will be elenoleel in a d.a .hoLir I' lC lormai, e.u. 12 1200, 
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however, that this upper-level wave is not associated with any well defined surface low, 
although a surface trough does extend from the Great Uikes to Kentucky with some reports 
of snow (Fig. 4). By 13/1800 (not shown), a second surface low has formed east of Cape 
Hatteras in response to the northeriy trough. 

The operational objective surface analyses without ERICA data do not accurately 
depict this developing cyclone situation. Thu.s, subjective analyses by Chalfant (1989), 
utilizing satellite imagery and ERICA aircraft and buoy data, are u.sed. At 14/0000 (Fig. 
6), three lows appear with central pre.ssures between 996 and 998 mb. Aloft (Fig. 7), the 
northern wave has moved off the East Coast of the US. Note that the isotherms are now 
m tre nearly in phase with the height ecmtcuirs. indicating little or no further baroclinic 
development aloft is expected. However, at the surface, the rapid intensification event 
begins at this time. 

The surface analysis at 14/0600 (not shown) shows the lows have merged and deep¬ 
ened approximately 17 mb over 6 h. By 141200 (Fig. 8). the low has plummeted to at 
least 968 mb (buoy at 37.5°N 63.5'\\'), a drop of 30 mb in 12 h. As w as expected from the 
14'0000 500 mb analysis of heights and tem|XTatures (Fig. 7). the height of the upper-level 
short wave (5340 m over Delaware) has slunvii little change tner the same 12 h period 
(Fig. 9). although a closed circulati('n has develi'petl. 

Operational numerical forecasts df the lOP-2 cyclone varied in accuracy. They arc- 
evaluated in detail by Chalfant (1980). None could successfully resolve the multiple low- 
centers at 14/0000. The NMC's Nested Grid McHlel (NGM) 48 h surface pressure forecast 
valid at 14/1200 (ncU shown) indicates cmly a 984 mb low, an error of 16 mb. However, 
the position forecast is only 100 km in error. The Global Spectral Model, also run at 
NMC, produced a much better forecast w ith central pressure of 964 mb (not show n) and a 
position error equivalent with the NGM. The operational NORAPS forecast (Simulation 1 
in section 5) correctly predicted the rapid deepening, slnwcing a 962 mb low at 48 h. but a 
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position error of approximately 350 km. However, as will be shown, the impressive fore¬ 
cast of rapid deepening at 48 h by NORAPS and the NMC spectral model is actually the 
result of poor forecasts at intermediate times. 
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IV. MODEL DESCRIPTION 


NORAPS is a regional model run operationally by FNOC, Monterey, CA over four 
areas: the Mediterranean Sea, the Indian Ocean, the North Atlantic Ocean, and the western 
North Pacific Ocean. The purpose of NORAPS is to provide high resolution forecasts of 
up to 48 h for areas important to Navy operations (Hodur 1987). The NORAPS model is 
easily relocatable to any area of the globe. The grid spacing is variable, both horizontally 
and vertical, but is limited by the speed and memory constraints of the computer system. 
For simulations in the current study, 60 km spacing is used on a 109 x 82 grid providing a 
model domain of approximately 6480 x 4860 km (Fig. 10). The vertical distribution of the 
21 vertical sigma levels is shown in Table 1. 


Table 1. NORAPS model sigma levels 


Mtxlel Level 

Siema 

Model Level 

Siema 

1 

0.010 

12 

0.600 

2 

O.O.T^ 

\? 

0.700 

3 

0.067 

14 

0.787 

4 

0.10.^ 

15 

0.850 

5 

0.150 

16 

0.890 

6 

0.200 

17 

0.920 

7 

0.250 

18 

0.950 

8 

0..”00 

19 

0.975 

9 

0..TS0 

20 

0.990 

10 

0.412 

21 

0.997 

11 

0.500 




NORAPS obtains initial data from radiosonde. PIBAL, tiircraft and ship observations. 
In certain regions, satellite vertical soundings and cloud-track w inds are used. Optimum 
Interpolation (01) is used for the analysis of these data with the previous 12 h forecast used 
;is the first guess field. In order t(' "spin up” the model, four 12 h update cycles are run 


10 







before a realistic simulation is performed. To begin the update cycle, the NOGAPS (Nava! 
Operational Global Atmospheric Prediction System) analyses are used as first guess field.s. 
Neither NORAPS nor the global model NOGAPS perform a moisture analysis. Thus, the 
previous 12 h moisture forecast is used as the analysis. 

Boundary conditions for the forecast fields in the NORAPS model run are obtained 
from NOGAPS forecasts. NORAPS uses one-way interactive lateral boundary conditions, 
which create a "blend zone" near the edge of the model grid. If the global model provides 
a poor representation of conditions on the NORAPS boundary, the forecast in the iruerior 
of the domain will eventually be contaminated. The length of time it takes for contamina¬ 
tion to occur depends largely on the speed at which the wind field advects the boundary 
conditions across the domain. The longer the forecast, or the smaller the model domain, 
the more likely it is that poor boundary conditions may adversely affect the NORAPS 
forecast. For the present study, the mtulel domain was situated so that its left boundary 
was far upstream of the incipient cyclone in carder to minimize the effects of the boundary 
conditions (Fig. 10). 

Topography in NORAPS is obtained from the Ntivy ten minute global data base. A 
terrain enveloping procedure developed by Wallace et al. (1983) is used to more effective¬ 
ly represent terrain effects in the model. A land-seti table with ten minute resolution is 
used by the model to define land and sea grid points. O'^ cr water, the albedo is set to 0.09, 
the ground wetness is 1.0, and the surface roughness is computed. Over ice. these values 
are 0.6, 0.0, and 0.0002435 m. respectively. Monthly climatokigictd values are used for 
the land grid points, as w'ell as the sea-surface temperature. 

The model uses an Arakawa Type-C staggered grid (Arakawa .and Ltimb 1977). Time 
integration is performed using a split-explicit scheme developed by Madala (1982). In this 
method, the solution is adjusted lor the linear terms which govern the fastest gravity 
modes, permitting use of a relatively large model time step. All :id\ection terms are 





second order. Fourth-order diffusion is included in the interior with second order diffusion 
on the first row and column of the grid. High frequency time oscillations are controlled by 
applying a Robert (1966) time filter with the smoothing coefficient set to 0.15. 

The parameterization of surface fluxes in the planetary boundary layer (PBL) are 
computed using the similarity theory (Louis et al. 1982), where matching solutions are 
obtained for the constant flux layer and the Ekman layer. The scheme is first order K- 
theory where K is a function of the bulk Richardson number. 

Since the lowest model sigma level above the surface is equal to 0.9975, a predictive 
equation for ground temperature given by Blackadar (1977) is used for land grid points. 
The radiation parameterization follows the methods of Sasamori (1968) for long wave 
radiation and Arakaw'a (1971) for short wave radiation. The effects of clouds are incorpo¬ 
rated into both radiation parameterizations, w’hile short wave radiation includes a diurnal 
cycle. 

Precipitation in NORAPS can occur from convective processes (sub-grid scale) or 
from large-scale condensation. Convection is parameterized using a Kuo-type scheme 
(Kuo 1965; Kuo 1974), in which precipitation is linked to low'-level moisture convergence. 
With each time step, the model determines whether a grid point meets the neces.sary condi¬ 
tions for convection. First, there must be moisture convergence in the PBL. Mathemati¬ 
cally, 

= [P.lg) jv q\dp > 0 

where is the moisture convergence, P^. is the surface pressure, g is gravity, q is the 
specific humidity, and the integral is performed over a depth of the atmosphere. 

Secondly, the value of 0^ above the LCL must be less than that at the LCL (i.e., a 
conditionally unstable layer must exist above the LCL). The depth of this conditional 
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instability must be greater than 300 mb. Finally, the relative humidity at the originating 
level must be greater than 80%. 


If the four conditions are met, the model transports heat and moisture vertically to 
each level (subscript k) using the equations given by Kuo (1965), 

^ _ Mt {Tc - T), _ 

' J/Cp(r.-r)rfp + ifLigc-g)dp 


^gk = 7 


(?c - g)k 


^ f Cp (Tc - T) dp -h j/L (gc - g) dp 


Here, T and q are temperature and specific humidity respectively, where the c subscript 
denotes cloud and unsubscripted denotes the environment. Cp represep*s specific heat and 
L is latent heat. Integrations are performed over the depth of the cloud (i.e., where (T^, - T) 
> 0 ). 

Kuo (1974) modified his parameterization with the inclusion of the "b” factor. This 
regulated the amount of heating verses moistening that occurs in the scheme. The fraction 
"b" is calculated followine 

100%-i? 

~ 100 % - 50 % ’ 

where R equals the integrated relative humidity over the depth of the cloud and 100% > R 
> 50%. Now the equations are modified as follows: 


(1 - b) Mt (Tc - T)k 
ifCp iTc-T)dp 


Sgk 


b Mt (gc - g)k 

j f ligc-g) dp 
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The role of "b” is to allow for rapid moistening initially when R is small (large b). As the 
column moistens, R increases and b decreases. This enhances the latent heating and de¬ 
creases the moistening rate. 

At the time of IOP-2 (December, 1^)88), the NORAPS model operationally used the 
Kuo (1965) sclicnie. In April of 1989, the Kuo (1974) scheme was incorporated. Thus, 
the simulations in section 5 reproduce the operational forecast of December, 1988, as well 
as the current operational version. 

Stable (or nonconvective) precipitation is assumed to occur when a grid point be¬ 
comes supersaturated. Moisture is condensed until saturation is once again reached. The 
condensed moisture is evaporated into the next layer below until it becomes saturated. 

The NORAPS mode! has been used operationally by the U.S. Navv since 1984. Skill 
scores show its forecast accuracy to be ctmiparable with other numerical models (Hodur 
1987). In addition, the NORAPS model has been used extensively as a research tool at the 
Naval Oceanographic and Atmospheric Research Ltihortitory and the Naval Postgraduate 
School. Recent studies (Liou et al. 1990; Wash et al. 1990) have utilized the model to 
study cyclogenesis during GALE. 




V. MODEL SIMULATIONS 


In order to investigate the physical processes of this rapidly deepening cyclone, 
numerical simulations of the cyclone were run using the NORAPS model. Forecasts of 48 
h were made beginning at 12/1200, 12 h before the initial development and 36 h before the 
rapid deepening. Since the intensification occurred with significant convection (Chalfant 
1989), the convective parameterization scheme of NORAPS was varied with all other 
physics in the model treated identically. The four convective parameterization schemes 
that were run are: 1) Kuo (1965) cumulus scheme (which was the operational version in 
Dpfpmhpr 1988): 2) Kuo (1974); 3) Kuo (1974) with vertical eddy diffusion; 4) convec¬ 
tive and stable precipitation turned off“. 

A. SIMULATION 1 

Simulation 1 of the IOP-2 event appears to begin well. A weak low forms at 3()°N off 
the Florida coast by 13/0000 (Fig. 11). But comparison with Fig. 3 shows that the forecast 
low is northeast of the verifying position by over 750 km. This is a large error for a 12 h 
forecast. 

Figures 12 and 13 present the (a) 1000 mb .streamlines tmd (b) 500 mb heights and 
vorticity for the analysis (12/1200) and 12 h forecast (13'0000) respectively. A weak 500 
mb disturbance is located over Florida in the analysis as indicated by the 8 s'^ vorticity 
contour (Fig. 12b). An examination of actual rawinsonde data (not shown) indicates that 
the disturbance did exist as a result ot cyclonic siiear across northern Florida.As this dis¬ 
turbance moves otf the coast, it becomes superimposed over the trailing cold front, as 


2. Hereafter, the Kuo (1%5) and Kuo (1074) schemes will be denoted KUOb.s tind KU()74 respectively. 




shown by the surface streamlines (Fig. 12a). This system develops rapidly during the 12 h 
period as the 500 mb absolute vorticity doubles from 8 s'^ to 16 s'^ (Fig. 13b). However, 
there has been no corresponding change in the height contours to account for such a 
change in vorticity. 

Note that the southerly short wave mentioned in section 4 is still over the southeastern 
U.S., as indicated by the vorticity maxima over Georgia in Fig. 13b. There is very little 
indication of the existence of this weak shortwave over Florida in the actual 500 mb analy¬ 
sis for 13/0000 (Fig. 2), although by 13/1200 (Fig. 5) it appears near 30°N 65°W, but 
apparently no '■■irface development took place. In other words, NORAPS incorrectly 
developed a surface low well to the east of the southerly wave in the IOP-2 event. The 
existence of this low at the surface is not verified by observational data, although the weak 
500 mb wave did, in fact, exist. Comparing Figs. 13;i and 13b, it is noted that the forecast 
1000 mb circulation is southwest of the vorticity maxima at 500 mb. Thus, the low is 
actually lying under negative vorticity advection. Closer inspection reveals that the low- 
level convergence zone northeast of the 1000 mb low is directly under the 500 mb vorticity 
maxima. 

Figure 14 shows the 12 h cumulative convective precipitation forecast by the model. 
Note that the precipitation maxima lies directly over the 1000 mb convergence zone north¬ 
east of the low (Fig. 13a). This is not .surprising, since l('w-le\ el mcristure ccmvergence is 
one of the four necessary conditions for convection in the Kuo scheme (see Section 3). 
The convective precipitation is also ccnncident with the 500 mb vorticity maxima (Fig. 
13b). While some surface obserxations dn indicate .some shower activity in this area, there 
is nothing to substantiate the predicted 2.8 cm of convective rainfall. Thus, the first 12 h 
of the model simulation can be summarized as follows; a weak upper-level short wave 
trough moved over a low-level baroclinic zone east of Florida, initiating cyclogenesis. 
Low-level convergence associateii with the surface low trieeeied convection northeast (d 






the incipient cyclone. The model convection produces an erroneous vorticity maxima at 
500 mb. 

This is even more apparent in Fig. 15 where the 18 h 500 mb heights and vorticity, 
and the 12 to 18 h convective precipitation fields are compared. The model continues to 
successfully predict the real short wave over the Georgia coastline (Fig. 15a). But it is 
obvious that the height and vorticity pattern to the east is anomalous. The model predicts 
12.1 cm of convective rainfall to occur at a grid point during the pa.st 6 h with a secondary 
maxima of approximately 7 cm to the east. Qmiparing this to Fig. 15a, note the colloca¬ 
tion of the 12.1 cm rainfall maxima and the small ridging in the 5640 m line at 500 mb 
likely due to the diabatic heating. Once again, the ccmvective precipitation lies over the 
convergent region north and east of the 1000 mb low- (not shown). 

As the upstream short wave in Fig. 15a merges with the developing system, rapid 
deepening takes place (7 mb in 6 h, 17 mb in 12 h. and 28 mb in 18 h). By 14/0000, the 
model has deepened the low to 980 mb (Fig. 16). A comparison w'ith the subjective analy¬ 
sis in Fig. 6 shows the failure of the model fc'recast. Only one low is predicted which is 
16-18 mb too deep and too far north, although a surface trough northwest of the low does 
hint at a secondary development as does the c(>nvective precipitation (not shown). The 
500 mb forecast (not shown) still contains the intense northerly short wave trough moving 
off the east coast of the US. How'e\ er. the initial low has already developed too much to 
allow' the new short wave to create a new low. 

Finally, Fig. 17 shows the 48 h NOF^APS forecast at the surface and 500 mb. The 
forecast surface low is only a few millibars too deep (assuming the 968 mb buoy is a true 
representation of the central pressure), but the position error is approximately 500 km (Fig. 
8). At 500 mb, the low- is similarly too deep ;ind too far n<'rth. 

In summary, it appear' that the (operational NOR.APS model predicted both 500 mb 
synoptic features reasonably w'ell; the southerly siK'rl w a\ e from Florida, and the norther- 
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ly wave from Gmada. Unfortunately, the model over-forecast ihe development of a weak 
low w'hich preceded the southerly short wave. By the time the northerly short wave ar¬ 
rived, the initial low' was too w'ell developed to permit secondary development. Instead, 
the northerly wave could only contribute to the deepening of the existing low. 

Several studies have show'n that enhanced precipitation occurs in an area extending 
from the northeast to the north and west of a surface low prior to the onset of rapid deepen¬ 
ing (Krishnamurti 1968; Johnson and Downey 1976; Kocin and Uccellini 1989). Kuo and 
Low-Nam (1990) found that this warm frontal precipitation w as crucial to the explosive 
development of the cyclone. The w arm front predicted by the model, as seen in the 1000 
mb wind field and the convective precipitation pattern, is in agreement with these findings 

Examination of the model output reveals that during the development of the surface 
low, the low is located either directly under the 500 mb vortieity maxima or upstream of it. 
One would expect the low to he loctited dow nstream t'f the vortieity maxima in the region 
of the positive v cirlicity advection. In additimi. convective precipitation charts indie;ite that 
the convective nature of this storm contributed heavily to its development. While ship 
reports and satellite imagery verify that convective activity was present (not show n), it is 
unlikely that the model simulation is accurate in this respect. To in\ estigate the conv ectiv e 
aspect of the model simulation, diagnostics of the eftects ig eon\ eetion are presented in the 
next section. 

B. SIMULATION 2 

A know'n deficiency of the KUOo5 scheme is that it added too much moisture aloft 
(Kuo 1974). In order to correct for this, the "b" factor was added w hich reduces the mois¬ 
tening w'hen near saturation is reached (see discussitm o1 consectiim pttrameieri/.ition in 
section .7). Simulation 2 used the KU074 .scheme to see if the change made any difterence 
in the lOP-2 c'se. 
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Figure 18 presents the 6 h cumulative convective precipitation predicted by the model 
valid at 13/0600. Comparing this with Fig. 15b, it is apparent that the KU074 scheme has 
reduced the convective rainfall amount by more than 505/. However, the surface and 500 
mb forecasts (not shown) at all times were nearly identical to that of Simulation 1(1-2 
mb differences in the central sea level pressure). The reason for the nearly identical fore¬ 
cast is that the stable precipitation was greater in Simulation 2. as shown by Table 2. In 
other words, since the Knn74 scheme caused less convective precipitation, there was 
more moisture available for the non-convecti\e scheme. As Table 2 indicates, the total 
amount of precipitativin i.n Simulation 2 is nearly equal that ('f Simulation 1. Therefore, the 
amount of latent heat released m the two simulations is nearly equ:il. resulting in compara¬ 
ble forecasts. Danard ( HhSs) and Luni and Fisberry (Ibis'/) trbserved a direct ct’rrehition 
between the maximum amount of convective precijutation and the pretiicted deepening 
r;ite. Kuo tmd l.ow-N'am (IbbO) huind that the amount of stable precipittition wtis key to 
determining the .amount of deepening. In the current stud> the tiual precipittition was 
ccirrehited to the deepening rtite. 

Table 2. Precipittition (cm) predicted by the NOR.-M’S model for Simulations 1. 2. 3. tiiid 
5. The amounts are ti 4 n h sum ot till grid points uith. precipitation ctuised b_\ the IOP-2 
cyclones. .Also included is the predicted centrtd pressure (mb) tit 4S 1,. 
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C. SIMULATION 3 


The inclusion of vertical eddy diffusion in the convection scheme was tested in 
Simulation 3. The effect of this diffusion is to mix the heat and moisture provided by the 
cumulus parameterization with the surrounding environment. This addition was expected 
to damp the effects of convection predicted by the mcHiel. As seen in Table 2, the convec¬ 
tive precipitation was further reduced by the eddy diffusion. How'ever, once again the 
stable precipitation was greater as a result. Thus. Simulation 3 showed little difference 
from the previous two simulations. 

D. SIMULATION 4 

.A fourth simuhition was conducted in which no precipittition v. ;is allowed. The goal 
of this simuhition was to ensure thtit the convective precipitation scheme was truly respon¬ 
sible for the pocir forecasts obserxed in the first three simulations. It was unrealistic since 
no con\'ecti\ e or stable precipitation was allmxed. In other words, when a grid point 
became super-saturtited, the model did nothing to correct this condition. Any development 
that takes place is due only to adiabtitiv •■'aroclinic dynamics, since there is no latent heat 
release. 

Figure F) shows the (;i) 1000 mb winds and (b) .'^00 mb heights :ind \orticity tor the 
12 h forecast. There tire mtirked differences between this simulatii'n and Simulation 1 
(Fig. 13). The initial short wave at 500 mb over f'lorida does not develop in Simulation 4 
(note that the enclosed area of <S s'^ vorticity netir 70 W in Fig. lOh is a minimum). 
Moreover, no devehipment from this system lakes place at KiOO mb (Fig. 10a). This 
indictites that even the small amount of convective precipitation (i.e. latent hetit release) 
generated in Simulation 1 (Fig. 14) is necessary for the development of the initial short 
wave. Figure Fhi shows the beginning id the low -lex el circulation along the Florida coast 
in response to the upper-level short w ave over (ieorgia. Thus, the mirdel has correctly 




ignored the initial w eak disturbance and developed a surface cyclone as a result of the 
second, stronger short w'ave. 

How-ever, the model forecasts in Simulation 4 are in error by not forecasting the 
vigorous development. The predicted surface cyclone is too shallow and too slow'. Figure 
20 shows a 1001 mb low at 14/0000 with an inverted trough extending along the US coast¬ 
line. Comparison with Fig. 6 demon.strales that p’-ecipitation was important, even at this 
early stage. The initial low' is 3-5 mb too shallow' and 475 km too far southwest. As one 
might expect, the rapid deepening of the surface cyclone is completely missed by the 
model in Simulation 4. The central pressure at 14/1200 is underforecast by 23 mb (not 
show'n). 

E. SU.MMARY OF SIMULATIONS 1 - 4 

Figure 21 summarizes the four numeucal simulations. It is seen that the model never 
exhibits a "rapidly deepening" phase, but rtither developed the cyclone at a constant rate, 
although predicted deepening rales betw een 13/0600 and 14/0600 (40 mb in 24 h) are 
certainly impressis e. In contrast, the observed cyclone exhibited the rtipid deepening in a 
6 h period'"’. The fact that, at 14/0000. the first three simulations were » 15 mb too deep 
and Simulation 4 w as only 5 mb kh' shallow, seems to indicate that the role of convection 
was rather small up to this point. Still, the errors early in Simulation 4 show that precipita¬ 
tion did play a role in the early stages of develcipment. Oh\ iously. though, the importance 
of latent heat release in the explosively deepening phase cannot be ignored. To investigate 
this more closelv, several diagnostic studies were performed on the model output and are 
discussed in section 6. 


'ITicrc is high confidence in the ;ictual ceniral pressure trace due to the tiirenifi and buoy me.'Lsurrrienis 
available. 





F. OTHER OPERATIONAL FORECASTS 


As mentioned in section 3, the NGM and Spectral models run at NMC produced very 
good operational forecasts in terms of position error (100 km at 48 h). However, the NGM 
was 16 mb too shallow in its central pressure forecast. The Spectral was much better (4 
mb too deep). But it too was similar to the NORAPS model in that it gave a very poor 36 
h forecast (Fig. 22). As with NORAPS, the Spectral model failed to resolve the multiple 
centers, and the central pressure of the low is 10 mb too deep. The 500 mb vorticity fore¬ 
cast is not as spurious as NORAPS, but since a smoothing procedure is applied to these 
operational charts, direct comparison is not possible. 

A comparison of precipitation forecasts from the NGM and NORAPS shows that the 
NGM was much more conservative in its prediction (not sht^wn). From 24 - 36 h, the 
NGM forecast a maxima of 1" (2.4 cm). while Simulation 1 predicted a maxima of 4.5" 
(10.9 cm) and Simulation 2 a maxima of 3.5” (8.5 cm). Precipitation forecasts for the 
Spectral model were not available. 





VI. MODEL DIAGNOSTICS 


In the previous section it was observed that w'hile the 48 h prediction by the NORAPS 
model was very good in terms of central pressure, the 36 h forecast (valid at 14/0000) was 
16 mb too deep. This was largely the result of intense convection early in the simulation. 
The model run with no precipitation showed only a small error in central pressure at 
14/0000. This section will investigate the effects of convection in the model simulations. 
In all figures of model data, the actual sigma level values were used to construct cross 
sections and Skew-T diagrams. 

A. TIME SECTIONS 

Section 3 described the Kuo .scheme of convection parameterization which is used in 
the NORAPS model. The effective result of convection is to transport latent heat and 
moisture from low levels (i.e., the PBL) upward. But before convection in the model 
begins, four criteria must be met (reviewcti in Section 3). 

It was already shown in section .s that the convective precipitation was occurring in 
areas with convergence at 1000 mb. implying thtit the criteria for low-level moisture con¬ 
vergence had been met. To see if the other three criteriti were correctly satisfied, time 
sections w'ere constructed at ERICA drifting buoy loctitions in which intense convective 
precipitation was falling in the model. Figures 23 - 25 show' the equivalent potential 
temperature, relative humidity, and mixing ratio time sections for Simulation 1 at 33°N 
7(EN (buoy 11369). The time sections clearly show that the Kuo scheme is accomplishing 
its goal. The profile of 0^ indicates a deep layer (surface to 850 mb) of conditional insta¬ 
bility due to cold air from the New England anticyclone flowing over the warm Gulf 
Stream w aters. The air is also nearly satiir.'ited up to 650 mb (Fig. 24). Thus, the neces- 




sary conditions for convection outlined in section 3 appear to exist. The potential tempera¬ 
ture time section (not shown) indicates weak static stability below 850 mb with stronger 
stability above this layer. Thus, the conditional instability (d0g/dz < 0) in Fig. 23 is large¬ 
ly due to the moisture profile. 

From charts of convective precipitation (see Figs. 14 and 15b), it was determined that 
very light precipitation from convection began in NORAPS between 12 and 18 h into the 
forecast simulation at the 33°N 70°W grid point. From hours 18 to 24, approximately 5 
cm of precipitation fell at this location. Figures 23 and 25 indicate that the convective 
parameterization began to take effect after about the 12 h mark, since the upper levels (700 
- 500 mb) begin to moisten creating a conditionally neutral layer aloft as shown by the 
vertical 0^ lines. The low-level conditional instability also begins to be removed after 12 
h. A sharp jump in the moisture profile occurs 18 h into the simulation (Fig. 25). the same 
time at which the significant convective precipitation begins. The entire column from the 
surface to 500 mb is nearly neutral by this time. The low-level moisture has been moved 
upwards as evidenced by the pocket of greater than 00% humidity in Fig. 24 and the drier 
70% contours below. After hour 21. the upper levels begin to dry rapidly (Fig. 24) indicat¬ 
ing the intrusion of drier air from the west, while the low levels begin to moisten as evi¬ 
denced by the 90% contour, setting up another conditionally unstable situation. The cause 
of this moistening is two-fold: evaporation of the precipitation into the drier low-level air 
and a shift in the low'-level wind direction from northeasterly to sttutherly (not showm) as 
the warm front passes by, advecting mtue moist air from the south. 

Figure 26 is a north-south cross section from a study by Reed and Albright (1986) of 
an explosively deepening cyclone in the eastern North Pacific. Note hinv the equivalent 
potential temperature contours are oriented vertically where convection is occurring. The 
0J. lines in the current study (Fig. 23) are also oriented vertically after 18 h as significant 
convective precipitation began to fall. In the cold air to the north (left) in Fig. 26. potential 


24 




instability is indicated by the d0g/dz < 0 up to 700 mb. This is similar to the left side of 
Fig. 23 where the cold air from the New England anticyclone is moving over the warm 
water. 

A comparison between the air temperatures measured by the buoy and the analysis 
temperatures at each model level was performed (not shown). This comparison revealed 
that the lowest model level (21 m) was 1.4°C warmer than the buoy. Gjnsidering that the 
New England anticyclone was advecting cold air over the warm water, one would expect 
to find colder temperatures above the surface layer. This anomaly was also observed at 
other buoy locations in the same region. The temperature difference became larger as the 
simulation progres.sed. The cause of this warming is mo.st likely the warm ocean waters. 
But a comparison of the measured sea-surface temperatures (SST's) and the model SST 
analysis showed good agreement. This appears Ki imply that the NORAPS model might 
possibly be modifying the cold air too rapidly. This warming by the warm ocean destabi¬ 
lizes the atmosphere, as seen in Fig. 23. 

Similar time sections were constructed for Simulations 2-4. It was pointed out in 
section 5 that Simulations 2 and 3 did not differ significantly from Simulation 1 in terms of 
the .sea level pressure forecast, but that the convective and stable precipitation predictions 
were different. The precipitation charts indicate that the total rainfall at the diagnostic grid 
point for Simulation 2 (KU074) was much le.ss than Simulation 1 (KU065), (« 1.5 cm 
compared with « 5.0 cm). The KU074 scheme intentionally decreases the moistening 
when saturation is nearly reached, which decreases the timount of precipitation. But as the 
moistening is decreased, the scheme c('rrespondingly increases the sensible heating. 
Additionally, stable precipitation is increased in the KU074 simulation. Thus, the overall 
heating profile is nearly identical to Simuhition 1 (KU065). 

This difference between the two convective parameterization schemes is shown in the 
time sections. Figures 27 and 28 are the equivalent potential temperature and mixing ratio 







time sections for Simulation 2. Comparing these to Figs. 23 and 25, both simulations look 
nearly identical for the first 15 h. But the conditional instability is removed more quickly 
in Simulation 1 as shown by comparing the vertical 0^ contours in Figs. 23 and 27. The 
moistening process is also slower in the second simulation, which is noted by comparing 
the 8 g/kg contour in Figs. 25 and 28. The first simulation (KU065) moistens the column 
more quickly than Simulation 2 (KU075), which removes the conditionally unstable layer 
faster and causes more precipitation. Thus, the intense convection in Simulation 2 is 
delayed as compared to Simulation 1. 

As would be expected, the time sections for Simulation 4 (no precipitation) are con¬ 
siderably different from the other three simulations. Figures 29 and 30 show the equiva¬ 
lent potential temperature and the mixing ratit' respectively. The initial iow-ltv’el condi¬ 
tional instability exists as in Simulations 1 and 2. However, with the convective parame¬ 
terization turned off, the model has no way to relieve this instability. Instead, the condi¬ 
tionally unstable layer deepens. Alstg note the marked differences in the moisture profiles 
of the twii simulations, with Simulation 4 indicating only ;i small amount of moistening in 
the low levels, and none aloft. During this period {15 to 27 h). the i(^w-level winds shift 
from northeasterly to easterly, again suggesting that ad\ection was the primary cause of 
the moistening observed in this simulation. 

B. IMPORTANCE OF INITIAL CONDITIONS 

From a qualitative analysis, the convective parameterization scheme appears to be 
performing correctly. The model atmosphere possesses the necessary crtnditions to cause 
convection and the Kuo scheme relieves the pi'tential instability by transpc'rting low -level 
latent heat and moisture upward. Httwever. the high values (tf convective precipitation 
noted in section 4 indicate that the model is perhrrming incorrectly in its prediction of 







convection. Since the predicted temperature and moisture play a large role in the cumulus 
parameterization, these fields are further investigated. 

The time sections in section 5a were coincident with a drifting buoy in the location of 
predicted convective precipitation. However, this did not allow for comparison with actual 
data at levels above the surface. In order to achieve this, a cross section was constructed 
using four Air Force dropsondes taken at 13/0000. Figure 31 is the equivalent potential 
temperature and mixing ratio cross sections along the line indicated on Fig. 14. The low- 
level potential instability created by the advection of cold air over warm water is clearly 
indicated by the "tongue" of low 0^ values near 900 mb. The corresponding NORAPS 12 
h forecast cross section for Simulation 1 is shown in Fig. 32. NORAPS also possesses a 
layer of potential instability. From Fig. 14 it is noted that convective precipitation has 
already been occurring in the region of the cross section. Thus, some of the instability has 
already been relieved, especially in the region between dropsonde 3 and 4. A similar cross 
section 6 h earlier (not shown) does in fact show that the model atmosphere was more 
unstable. 

Comparison of Figs. 31 and 32 reveals a disturbing fact: the model atmosphere is 
much moister in the low levels, typically by 2 g/Tcg. Analysis and 6 h cross sections (not 
shown) show little change in the moisture field. Thus, it appears NORAPS was initialized 
with a poor moisture field. As was previously mentioned. NORAPS does not perform an 
actual moisture analysis. Instead, it uses its previous 12 h moisture forecast as the analy¬ 
sis. Obviously, this could result in a situation such as that in the present study. 

Figure 33 shows a comparison of the actual Air Force dropsonde 3 and the corre¬ 
sponding NORAPS 12 h forecast sounding (denoted by /-\F03 in Figs. 31 and 32). While 
the overall soundings look similar, closer inspection reveals that NORAPS is both moister 
and warmer than the observed atmosphere. Also note that the observed inversion between 
850 - 700 mb (Fig. 33a) is not predicteti by the model (Fig. 33b). The result is a more 
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unstable sounding in the model, as indicated by the stability indices. Similar soundings for 
all other (18) dropsondes show that the model is more unstable than the real atmosphere. 
Tills is somewhat surprising since the actual atmosphere "looked" more unstable than the 
model in the 0g cross sections (Pigs. 31a and 32a). ALso, convective activity in this region 
had already served to stabilize the model atmosphere somewhat. Clearly, the model is not 
doing a good job of predicting the temperature and moisture fields in the lower levels. 
And since the analysis fields (not shown) are even more unstable, it is suspected that the 
initial NORAPS fields poorly represented the actual atmosphere. 

In order to determine if the poor representation of the temperature and moi.sture was 
due the lack of data, comparisons were made between the model and the 12/1200 sound¬ 
ings (not shown) at West Palm Beach, FL (72203). Charleston, SC (72208), Cape Haiteras. 
NC (72304), Nassau, Bahamas (78073). and Bermuda Island (78016). In all cases but 
Bermuda, the model analysis was nearly identical to the observaition data. As with the Air 
Force dropsondes. NORAPS w'as too vxiim and moist ;it Bermuda. Therefore it appears 
that the NORAPS analysis gives a correct depiction of the titmosphere when data are avail¬ 
able, but may be poor over the data sparse ocean regions. Model 12 h forecast errors are 
passed on to the next analysis through the first guess field if there is no data to refute it. 
Thus, a poor 12 h forecast can result in a pentr analysis. 

It was mentioned in section 3 that for the present study, four 12 h forecasts had been 
made in order to simulate the update eycie m an ('peratii'nal ens ironment. Hodur (1987) 
showed that the inclusion of a 12 h update cycle improved NORAPS forecasts. Thus, the 
NORAPS analysis could differ from the gh'bal mcnlel (NOO.APS) analysis because of the 
use of the previous 12 h forecast :is the first guess field. Also in section 3 it was pointed 
out that neither NORAPS nor NOOAPS analyze moisture. Rather, the initial moisture 
fie.d for a model simulation is the previous 12 h torecast. Surface tempertitures measured 
by ships and buens are also not used in the analyMs. 





To investigate the impact of these procedures. Fig. 34 compares the 2 m air tempera¬ 
ture analysis from NORAPS (34a) and NOGAPS (34b) with actual ship and buoy observa¬ 
tions. The ERICA buoy data shown in this figure was not used in the model analysis. 
Note that south of 35°N, NORAPS (Fig. 34a) is generally warmer than the observations, 
by as much as 6°C. The NOGAPS analysis (Fig. 34b) shows much less detail in the 
temperature field, especially near the coastline. However, the analysis appears to be in 
much better agreement with the observed temperatures, although the ERICA buoys near 
32°N are still colder than the analysis. Since the temperature is not analyzed directly by 
either model, these differences are largely a result of the previous 12 h forecast. 

Figure 35 is the initial 1000 mb vapor pressure (moisture) fields for the NORAPS 
(35a) and NOGAPS (35h) models. The NORAPS model again shows more detail in the 
1000 mb moisture field than does the global model. However, the regional model is 
considerably moister, especially south of 35°N (note the kK’ation of the 15 mb contour as it 
crosses 70°W). Comparison of actual moisture values at 12/1200 (not shown) to the 
NORAPS and NOG.APS analyses show that NOGAPS is more much more realistic, and 
may in fact be slightly too dry. The better moisture analyses in NOG.APS are actually 
better 12 h forecasts, since there is no actual analysis of moisture. 

The causes for the differences between the NORAPS ;ind NOGAPS low-level temper¬ 
ature and moisture fields are beyond the scope of this study. As mentioned in section 6a. 
one possibility is that as the continental pohir air mt'ves over the Gulf Stream waters, the 
NORAPS model is modifying the air too rapidiv, thus mtiking it warmer and moister in the 
low levels. The drier NOGAPS fields might also be the result of different model physics, 
such as the use of the Arakawa-Sehiibert (1M74) scheme bm convection (Rosmond'^. 

4. Dr. Thomas H. Rosmond is the head of the IVedietiim Systems bnineh of the Ntival Oeetinogniphic & 
Atmospheric Research I.abor;itory, Atmospheric Directorate. Monterey. ('A. 




personal communication), rather than the Kuo (1974) scheme used in the NORAPS model 
for the current study. 

C. SIMULATION WITH NOGAPS INITIAL FIELDS 

In order to test the hypothesis that the excess temperature and moisture in the 
NORAPS initial fields are adversely effecting the model simulations, a fifth simulation 
was carried out where the model was initialized with the global analysis fields (i.e., no 
regional update cycle). The model for the this simulation used the KU074 convection 
scheme. Additionally, minor improvements were made to the boundary conditions, but 
these have relatively small influence on the forecast (Liou, personal communication). 

A subjective comparison of the NOGAPS surface and 500 mb analyses with the fields 
used to start the first four simulations showed scmie minor differences, but these were 
primarily out of the region of interest. The important weak disturbance over Florida was 
depicted in both analyses by the 8 s'* vcuticity contour. Thus, the major difference be¬ 
tween this simulation and Simulation 2 is the ini'ial low-level temperature and moisture 
fields. 

The cumulative convective precipitation from 12 to 18 h is shown in Fig. 36. Com¬ 
paring this with Fig. 19 from Simulation 2. it is seen that the drier initial moisture iicid has 
had an effect on the convective precipitatii'n. both in amount and distribution. The corre¬ 
sponding sea level pressure fields (not shown) differ only by ] mb. Likewise, the 500 mb 
fields show little difference, although the new' simuhition moves the .short wave somew'hat 
faster. Still, the model is developing the w rong inititil low. But by hour 30 (Fig. 37), an 
improvement is noted. Since the inititil low at this time is not as deep as in Simulation 2 
(998 mb vs 991 mb in Simulation 2, not show n). NORAPS attempts to develop the sec¬ 
ondary low to the northwest. The subjective analysis from Chalfant (1989) at this time 
(not shown) shows the ob.served second low in gocnl agreement with the forecast position 
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and central pressure. However, the primary low has a large error in both position and 
central pressure, although it is closer to the ob.served 1002 mb value than Simulation 2. 

By 14/0000, the two simulations are considerably different. Figure 38 is the surface 
pressure forecast for Simulation 5. (The forecast by Simulation 2 at this time is identical 
to Simulation 1 (Fig. 16) and is thus not shown). The 989 mb low is still in error by 7-9 
mb. But this is a considerable improvement over Simulations 1-3. Note that the low in 
Simulation 5 is moving much more rapidly. Anthes et al. (1983) point out that latent 
heating reduces the cyclone's speed of movement. As Table 2 indicates, this simulation 
produced much less precipitation, thus agreeing with this finding. 

Table 2 also reveals the effect of the reduced precipitation. The 48 h central pressure 
forecast is 971 mb, in close agreement with the obsers ed 968 mb. Note tiiat it appears that 
the amount of predicted deepening is related to the loltil amount of precipitation, not the 
convective or stable amounts alone. KU074 reduced the amount of convection from that 
of KU065, but the stable precipitation incretised resulting in :i similar amount of latent 
heating. Only when both types of precipitation were reduced was the deepening rate 
improved (see Fig. 21). 

Therefore, it appears that the initial h'w -level temperature and moisture fields play a 
significant role in the de\ elopment of the cyclone. The cumulus parameterization is sensi¬ 
tive to small changes in the moisture field ;md correspondingly varies the amount of pre¬ 
cipitation, and thus the timount of latent heat release. The result is a difference in both the 
track and intensity of the predictetl cyclone, w ith both improved over the earlier simula¬ 
tions w'ith the NORAPS update cycle. The ex.act ctiuse of the excessive heat and moisture 
in the NORAPS analyses is beyond the sco|ic c>f the current study, but is most likely linked 
to the boundary layer parameterization (mixing) and vertical grid resolution (see Table 1). 
Again it should be noted that this is just one ctise :ind th:it similar studies mtty show results 
which support or contradict these findings. 



VII. SUMMARY AND CONCLUSIONS 


Considerable work has been done in the area of numerical prediction of rapid oceanic 
cyclogenesis. Several studies have focu.sed on the latent heating aspect as the key mecha¬ 
nism controlling the deepening rate. Yet, in nearly all of these studies, the deepening rate 
was typically less than observed. While the operational NORAPS 48 h forecast valid at 
14/1200 UTC appeared quite good, a closer in.spection reveals that the intermediate fore¬ 
casts were rather poor, with NORAPS over-forecasting the development of the ERICA 
lOP-2 cyclone. 

An investigation into this numerical forecast showed that the NORAPS model incor¬ 
rectly developed a low too early in the forecast due to pr(4'»lems in the convective parame¬ 
terization. Furthermore, the development of this low was too rapid, with an over-forecast 
central pressure of 18 mb at 14/0000 UTC. Unrealistically large amounts of convective 
precipitation appear to be responsible for the overdevelopment, as the convection generat¬ 
ed spurious amounts of vorticity at 500 mb with no appreciable change in the height pat¬ 
tern. 

Varying the cumulus parameterization scheme decreased the amount of convective 
precipitation, but since the stable precipitation c('rrcspimdingly increased, tb.is resulted in 
very little change in the surface pressure :md 500 mb height forecasts. Danard (1985) 
found a direct correlation between the mtiximum amount of convective precipitation and 
the predicted deepening. This fact was not observed in the present study. Kuo and Low- 
Nam (1990) emphasized the importance <4'warm frontal sttible precipitation to the rate of 
deepening. They found that the simulations using the Arakawa-Schubert convective 
parameterization were superior to those using Kuo. They concluded that this was primari¬ 
ly the result of the former scheme precipitating much less in convection than the latter. 





leaving more available moisture for the stable precipitation. Stable precipitation is typical¬ 
ly more efficient than convective precipitation at deepening a cyclone since it occurs lower 
in the atmosphere. As Anthes et al. (1983) points out, the lower the level of maximum 
heating, the more intense the resulting surface cyclone. For the current study, the total 
amount of precipitation was the only link to the predicted deepening rate. However, the 
importance of the location of the precipitation, both horizontal and vertical, cannot be 
overlooked. 

In order to verify that the latent heat release was the cause of the overforecast deepen¬ 
ing, a simulation with no precipitation was run. This alteraticm dramatically slowed the 
development of the initial cyclone, but failed to adeqiiatelv predict the explosive deepening 
w'hich was observed. The underprediction of the low prior t(^ the rapid deepening phase 
indicates that precipitation w’as important in the earlv sttiges. 

Detailed diagnostics were performed t(' tmalyze the role of the convective precipita¬ 
tion. It w'as determined that the cumulus ptirameteriztition scheme was working properly. 
How'ever. comparisons of model analyses .and forectists to ob.served data indicate that the 
model was poorly initialized, especitilly in low levels tner the ocetin. resulting in a great 
deal more potential instability. The comptirisims w ith HRICA dropsimdes also showed 
that the model forecasts w ere too w arm in the low levels, w hich mav be a result of a poor 
surface flux parameterization. It w as reci'gnized that the lack of a moisture anaUsis also 
might be partly responsible for the poor initialization. In ('rder to investigate this, a fifth 
simulation was conducted which used the global analvses as inititil fields. These global 
fields were found tir be drier anti more realistic. The resulting forecasts were improved, 
both in sea level pressure and precipitation amounts, alttanigh the model still deveh^ped the 
wrong low', and then failed to deepen it enough. 

A comparison with tw'o NMC models, the Nested Cirid Model and the Global Spectral 
Model, showed simihir deficiencies t<' those fountl in NORAPS. Both mtulels developed 




the initial low and tailed to resolve the secondary development off the Virginia - North 
Carolina coastline. While the NGM model was l(i mb too shallow at 14; 1200, the spectral 
model actually overforecast the development by 4 mb, similar to the operational NORAPS 
forecast. Both models did very well in predicting the 48 h location of the cyclone. 

This study agrees with previous findings on the importance of latent heating in the 
development of wintertime oceanic extralropical cyclones. The quality of numerical 
forecasts relies on the quality of the initial data. In data sparse regions, model forecast 
errors often propagate through successive forecasts via the first guess field unless there are 
available data to correct the analysis. Parameterizations f(U convection and surface fluxes 
are very sensitive to details in the miHlel atmosphere. Incorrect specification of this initial 
atmosphere can have an ad\ erse effect on the forecasts even if the parameterization 
schemes are correct. Thus, the common lack (d'data is :m important problem in the correct 
prediction of rapidly developing cycUmes. 

Since the initial moisture field tippears to have had a large effect on the NOR.APS 
forecasts, future work should be done in this ;irea. The incorpirration of an actual moisture 
analysis, rather than the pres ii'us 12 h forecast, cmikl have a significant impact on explo¬ 
sive cyclogenesis forecasting. Additionally, it was noted that the low -level temperatures in 
NORAPS were often too warm in data sparse regions, further adriing tir the potential insta¬ 
bility in the model. The current studv used onl) tiu'se I'bsersations which were ('per.ition- 
ally available at the time (4 lOP-2. The extra buox, rawmsonde. and dropsonde data col¬ 
lected as a part of ERICA should be mtegr.iterl into the NOR.APS analyses to iinestigaie 
the importance c.f these data in numerical forecasting. 
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Fig. 1. FNOC 500 mb analysis for 1200 UTC 12 December 19S8. Geopotential height 
(solid) IS in dekameters (+64 represents 564 dam) and temperature (dashed) is in '"C. 















f-'ig. 3. f'NOC surtacc analysis for (K)OO ll'I'C' 13 December 19NiS. Sea le\el pressure 
contours are in mb minus lOOO mb ( + 20 represents 1020 mb). Surface re|)orl iiuls aie 
in kt and temperatures are in °C. See text for explanation of circled buoy. 
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in Fig. 1 except for 1200 UTC 13 December 1988. 
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Fig. 6. Subjective surface analysis for 0000 UTC 14 December 1988. Surface pres.sure 
contour interval is 2 mb. Surf ace observations are of sea level pressure (125 = 1012.5 mb), 
temperature (in °C) and winds (in kt). From Chalfant (1989). 
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Fig. 8. As in Fig. 6 except for 1200 UTC 14 December 1988. 

















Fig. 9. As in Fig. 1 except for 1200 UTC 14 December 1988. 
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Fig. 11. 12 h forecast of sea level pressure (contour interval 4 mb) and 
1000-500 mb thickness (dashed, contour interval 60 m) for Simulation 1. 
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Fig. 12a. Initial 1000 mb streamlines for Simulation 1. 
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Fig. 12b. Initial 500 mb geopotential heights (contour interval 60 m) and 
absolute vorticity (dashed lines, contour interval 2 s'^) for Simulation 1. 






Fig. 13b. As in Fig. 12b except tor 12 h torecast. 
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Fig. 14. 12 h forecast of cumulative 12 h convective precipitation (contour 
interval 1 cm). Solid line indicates cross section location in Figs. 31 & 32. 
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Fig. 15a. 18 h forecast of 500 mb geopoteniial heights (contour interval 60 
m) and absolute vorticity (dashed fines, contour interval 2 s'^) for Simula- 






Fig. 15b. 18 h forecast of cumulative 6 h convective precipitation (contour 
interval 1 cm) for Simulation 1. 
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Fig. 16. As in Fig. 11 except for 36 h forecast. 
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Fig. 17b. As in Fig. 12b except for 48 h forecast for Simulation 1. 
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Fig. 18. As in Fig. 15b except for Simulation 2. 
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Fig. 19a. As in Fig. 13a except for Simulation 4. 
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Fig. 19b. As in Fig. 13b except for Simulation 
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Fig. 21. Observed and f jrecast central pressures (CP) in mb. Observed CP taken trom 
Chalfant (1989). Forecasts are for Simulations 1-5 (see text for explanation on simuki- 
tions). 
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Fig. 22a. NMC Global Spectral Model 36 h forecast of sea level pressure 
(contour interval 4 mb) and 1000-500 mb thickness (dashed lines, contour 
interval 60 m). 
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Fig. 22b. NMC Global Spectral Model 36 h forecast of 500 heighi.'< (con¬ 
tour interval 60 m) and vorlicity (dashed lines, contour interval 2* ). 


62 









YrfTi ["I 



























Fig. 28. As in Fig. 25 except for Simulation 






























Fig. 31. Vertical cross section of (a) equivalent potential temperature (contour interval 2 
K) and (b) mixing ratio (contour interval 1 g/kg) of four Air Force dropsondes at UtKK) 
UTC 13 December 1988. See Fig. 12 for location of cross section. 
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Fig. 32. As in Fig. 31 except for Simulation 1 12 h forecast. 
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Fig. 33. Skew-T Log P diagram of (a) Air Force 03 dropsonde and (b) corre¬ 
sponding NORAPS 12 h forecast valid at 0000 UTC 13 December 1988. 
Temperatures are in °C. Stability indices at right are Showalier Index (SI). 
Lifted Index (LI), K Index (KI) and Totals Index (TTl). 












Fig. 34a. 1200 UTC 12 December 1988 analysis of 2 m air temperature by 
NORAPS. Differences between observed ship and buoy temperatures and 
the model analyses are also plotted. A positive value indicates that the 
analysis is w'armer than the observation. Underlined numbers indicate 
ERICA buoys. 









Fig. 34b. As in Fig. 34a except for NOGAPS analysis. 













Fig. 35a. 1200 UTC 12 December 1988 analysis of 1000 mb vapor pres¬ 
sure (contour interval 5 mb) by NORAPS. 
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